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Abstract 



We calculate the minimal surface bounded by four-sided figures whose projection on 
a plane is a rectangle, starting with the bilinear interpolation and using, for smooth- 
ness, the Chebyshev polynomial expansion in our discretized numerical algorithm 
to get closer to satisfying the zero mean curvature condition. We report values for 
both the bilinear and improved areas, suggesting a quantitative evaluation of the 
bilinear interpolation. An analytical expression of the Schwarz minimal surface with 
polygonal boundaries and its 3-dimensional plot is also given. 
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1 Introduction 

In mathematical modeling it is not uncommon to need a surface that spans a 
known boundary and has the least value of a related quantity, say, area. If the 
least area is desired, the problem is termed in the mathematical literature as 
the Plateau problem, namely minimizing the area functional 
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Here Q G is a domain over which the surface X is defined as a map, with 
the boundary condition X{dfl) = F. It is known [1] that the first variation of 
A{X) vanishes if and only if the mean curvature if of X is zero everywhere in 
it. Thus to get a minimal (or, more precisely, a stationary) surface, we have 
to solve the differential equation obtained by setting the mean curvature H 
equal to zero for each value of the two parameters, say, u and v parameterizing 
a surface spanning the fixed boundary. In a numerical work, the problem 
has to be discretized by choosing a selection of the numerical values of the 
two parameters and finding the minimal- surface-position for each pair of 
the values. If the given boundary is a four-sided figure whose projection on a 
plane is a rectangle, the surface positions become simply the heights above the 
Mu-plane. Such 'numerical heights' and resulting 'numerical minimal surfaces' 
have been computed in ref. [8] for a variety of closed curve boundaries. 

In this paper, we report, in the section H] below, a modification to their algo- 
rithm that uses linear combinations of the Chebyshev polynomials as heights 
at the discretized wf-positions. In this way, we replace in the algorithm arbi- 
trary heights by linear combinations of convenient polynomials with arbitrary 
coefficients. 

The immediate advantage of this use of polynomials has been a reduction in 
the discretization error and a better convergence. Polynomials are (smooth) 
analytic functions having simply calculable derivatives. We have also carried 
further our efforts to find analytic surfaces that can be taken as 'approximate 
minimal surfaces': 

1) We read initial heights from a ruled analytic surface spanning our fixed 
boundary, namely the bilinear interpolation introduced in the section [31 And 
for knowing how much heights changed through our numerical minimization 

2) we compared the areas of the numerically found points (or 'numerical min- 
imal surface' explained in section 14.21 with those of the bilinear interpolation 
for each of the selected boundaries. 

Through this quantitative comparison, something missing in the previous 
works, we suggest to a user of a minimal surface bounded by four straight 
lines a prescription that may well save almost all the computer programming 
and CPU time spent in implementation, say, the algorithm of refs. [8]: the 
approximate equality of the areas of the bilinear interpolation and numerical 
minimal surface strongly suggests that the simple bilinear interpolation itself 
may work as a 'minimal surface' for many mathematical models that need 
minimal surfaces bounded by four straight lines. 

The only ruled surface, other than the plane, which is a minimal surface is 
a helicoid [Ij. As one boundary of a helicoid must be part of a helix, which 
is not a straight line, the boundary of a helicoid cannot be composed of four 
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straight lines. In this way there cannot be at least a ruled surface which is a 
minimal surface bounded by four straight lines. 

Since the calculation of the area given by the surface coordinates would be 
possible only numerically, it is technically important how to evaluate the min- 
imal surface area accurately, and evaluate deviation from the ruled surface 
whose area can be evaluated analytically. An area of bilinear interpolation is 
to be compared only with the numerically calculated 'minimal surfaces'. (See 
the section 14.21 below for a description of the algorithms we used to calculate 
areas of the 'numerical minimal surfaces' along with the resulting numerical 
values.) 



2 Plateau problem 



For a locally parameterized surface X = 'X.{x,y, z{x,y)), the mean curvature 
H is defined as 



rj 9llh22 - ^9l2hl2 + 922hll 

^ = 2 ' 

911922 — 9l2 



where 



9ii = {Xu,Xu), gi2 = {Xu,X^) and g22 = {Xi„X^) 

(3) 

are the 1st fundamental form and 

hu = {N,Xuu), hi2 = {N,Xuv) and /i22 = {N,X^^) (4) 

are the second fundamental form. Here 



iV= l^""^"l (5) 
is the unit normal of the surface. 

The vanishing condition of the numerator of H becomes 

F( "^^^ ll ( 2^^*^^ ll f'^'^Vl f6) 

dy'^ \ \dx I j dxdy dxdy dx'^ \ \dy I 
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We are interested in evaluating the area bounded by skew quadrilateral [TO] 
whose boundary is composed of four non-planar straight lines connecting four 
corners Xqo, Xqi, Xio and Xn. 

The Plateau problem for polygonal boundaries was studied by Schwarz, Weier- 
strass and Riemann [5|7P^ . 

The minimal surface whose bounding contour is the skew quadrilateral consist- 
ing of four edges 0, 5(0, -i, -^), C(-i, 0, ^) and D{0, i, 
was calculated by Schwarz [7] using the Weierstrass-Enneper representation. 
An extensive derivation of the minimal surface is given in [SP] . 

In this theory, every simply connected, open minimal surface with normal 
domain 11 is shown to be expressed in the form 

r = r{a,(3) = ro + Re r F{'y)d'y;'y CU (7) 

Jo 

where -^(7) is a non- vanishing analytic vector in 11 satisfying F'^ = (pK'j) + 
01(7) + 01(7) = 



One works with $(7) = ^ (Ml) -iMl))/'^ and ^(7) = ^(Ml) + ^Ml))/'^ 
and 2<I>\E' = (/)3 

When $ and do not have the common zero, the following expression was 
obtained: 

x = xo + Re (<l>^ - ^^)d7 

JO 

z = zo + Re r 2md-f (8) 
^0 

Using the mapping ^{■j) = \l'(7)/$(7), and defining $(7)^(i7 = R{u)du!, 
Schwarz obtained the expression 

x = Re J'^{l-uj'^)R{uj)duj 

/CO 
{l + iu^)R{uj)duj 

z = Re r 2uR{uj)duj (9) 



where 

R[u) 
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The integral can be done analytically, whose detail is given in the Appendix. 



3 The Bilinear Interpolation: 



We try to approach the minimal surface for the boundary composed of four 
non-planar straight lines connecting four corners xqo, xqi, Xiq and Xn by im- 
proving upon a surface that spans this boundary, namely a hyperbolic paraboloid 

m 



x(m, v) = [1 — u u] 


Xqo Xqi 




l-v 










XlO Xii 




V 



(10) 



(Hyperbolic paraboloid is a bilinear interpolation; it might interest the reader 
that this is a special case of the general bilinear interpolation, termed the 
Coons Patch For the corners we chose, for a selection of integer values of 
d and r: 



xoo = ri 



xio = r4 



xoi = rs 



xii = r2 



We consider two types of configurations of the four corners: ruledi and ruled2. 
In the case of ruledi we choose 

ri = (0,0,0) r2 = (r,rf,0) rg = (0, rf, rf) r4 = (r,0,rf). (12) 

The mapping from (u, v) to (a;, z) in this case is 



x{u, v) = ur 
y{u, v) = vr 

z{u, v) =ud + vd{l — 2u) 



(13) 



In the case of ruled2 we choose 

ri = (0,0,0) r2 = (r,r,0) r3 = (0,r,c/) 
The mapping from {u, v) to (x, y, z) in this case is 



r4 = (r,0,rf). (14) 



x(m, v) = ur 
y{u, v) = vd 

z{u, v)=ud + vd{l — 2u) 



(15) 
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Fig. 1. The ruledi surface (r=l,d=2). The horizontal plane is expanded by y, z and 
the height is x. 




Fig. 2. The ruled2 surface (r=2,d=l). The horizontal plane is expanded by x,y and 
the height is z. 

These definitions are such that for r = d the four position vectors lie at the 
corners of a regular tetrahedron. The Fig. 1 and Fig. 2 below are 3D graphs 
of the hyperbolic paraboloid for a choice of corners mentioned in eqs.f ll2p and 

m- 

For a surface to be minimal, its mean curvature vanishes everywhere [1]. The 
expression for the mean curvature, calculated using eq.([2]) of our bilinear in- 
terpolation is 



-2dr{l-2u){l-2v) 



for the ruledi and 



[r2 + r2(l - 2m)2 + d\l - 2t;)2]3/2 ' 



-2r3(l - 2m) (1 - 2v) 
d[d^ + r2(l - 2m)2 + r2(l - 2t;)2]3/2 ' 



(16) 



(17) 



for the ruled2. 
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Fig. 3. The numerically fitted ruledi surface (r=l,d=2). 

The mean curvature for the surface is zero only for the u = ^ line and the 
f = I line, whereas for a minimal surface this should be zero for all values of 
u and V. 

4 The Numerical Work: 

The solution of the Plateau problem was formulated by Courant[2] as mini- 
mization of the Dirichlet integral 

Ed{u) = \ I |Vn|2 
2 Jn 

where |VmP = iii^dudu)^ where du is the matrix of partial derivatives of u in 
an orthonormal basis. In [13j, a mapping to the conjugate minimal surface was 
considered in the minimization process. In [9], a diffeomorphism Mq : — > R'^, 
where VL C R^, UQiVt) C 5(0) and maps u{-,t) : Vt R^, (t > 0), which 
satisfies 

— -Asit)U = infix (0,T) 
with appropriate Dirichlet boundary condition was considered. 

In [8], more direct minimization of the numerator of the mean curvature H 
using parallel computer was performed. In the generalized Newton's method, 
the minimization of F{z) of eq.Q is achieved by the iteration 

^k+i ^^k_ DF{z'')~^[F{z% (18) 

where DF{z^)~^ is the inverse of the functional derivative that satisfies 

DF{z^) [z''+^ - z^] = -F{z^). (19) 
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We consider (A^ + 1) x (A^ + 1) lattice grid points {ui, vj), {0 < i < N,0 < j < 
N) and corresponding z{ui, Vj). We keep same number of grid points indepen- 
dent of r and d. In the discretized system z'^~^^{ui, Vj) is defined from z^{ui^ Vj) 
by adding dz^~^^{ui, Vj) which can be calculated by solving the linear equation 
expressed by a matrix C defined by the first and the second fundamental form 
as 

Cadz''+\u, v) = -Fiz^u, v)) (20) 

Businger et. al. [H] gave a Mathematica code to define the matrix C. In our 
problem of improving the surface starting from the bilinear area, the dis- 
cretization error in the replacement like 

dz_ ^ z{ui+i,Vj) - z{ui,Vj) 
du du 

is large and the convergence was poor. 

The reason would be lack of explicit third order polynomial term in the evalu- 
ation of dz^^^ in the numerical methods which manifests itself in the fact that 
C(j_ij) and C(j+ij) are identical. Thus we evaluate the first and second fun- 
damental form on the discretized system by using the Chebyshev polynomial 
expansion [B]. 



4-1 Chebyshev Polynomial Expansion 

The Chebyshev polynomial of degree n is denoted Tn{x) and is given by 

Tn{x) = cos(ncos~^x) (22) 

where the range of a; is [—1, 1] and their explicit expressions are given by the 
recursion 

To(x) = l, T„+i(x) = 2xr„(x)-T„_i(x) (n > 1) (23) 
The zeros of T„(x) are located at 

Xfc = cos(^^^^^^), {k = 0,l,--- ,n-l) (24) 
n 

If Xk {k = 0,1,-- - ,m— 1) are the m zeros of Tm{x), the Chebyshev polynomial 
satisfies the discrete orthogonality relation for i,j < m, 
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m— 1 

J2 Ti{xk)Tj{xk) = 

fc=0 



i^j 

m/2 i^jy^O (25) 
m i — j — 



We first map z{ui,Vj),{0 < Ui < 1,0 < vj < 1) onto z{xi,yj), {—1 < Xi < 
1,-1 < Hj < 1) and interpolate values at zeros of the Tisf+i{x) defined as xi 
(/ = 0, 1, • • • , A^) and TN+i{y) defined as i/m, (m = 0, 1, • • • ,N), i.e. z{xi, ym)- 

We define c{xi, n) {n — 0, • ■ • , N) as 

2 

c{Xh n) = 51 z{Xhym)Tn{ym) (26) 

m=l 

and interpolate at y = yj via 

^(^^ %) = ^(^'' ^)Tn{yj) - lc{xi, 0) (27) 
Partial derivative in y is performed by replacing T„{y) by 



dy 



Tn{y)9yz{xhyj) = ^c(a;;,n)T^(%). (28) 



So far the x-coordinate is restricted to zero points x''. Now, interpolation to 
X = Xi is performed by 



2 

yj) = ^(^'' yj)Tn{xi) (29) 



We define also dyc{n,yj) as 



z=i 



dyc{n,yj) = j^^-p^ 9yz{xi,yj)T^{xi) 
The values on the mesh points z{xi,yj) are 



(30) 



yj) = Yl yj)^n(2;i) - ^c(o, %■) (31) 

and the derivatives dxz{xi,yj) and d'^z{xi,yj) are 

^x5(a;i, %■) = Y yj)Tn{xi) (32) 
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dlz{xi, Vj) = J2 c{n, yj)T'^{xi) (33) 

n 

dxdyz{xi, Vj) = dyc{n, yj)T^{xi) (34) 

n 

In the linear equation 

Cdz'+\ij^ = (35) 

the matrix C in the left-hand side(lhs) is a sparse matrix that contains at least 
nine non-vanishing elements in each row. Around the position (0 < i < 
N,0 < j < N) the elements for the nine nearest neighbors of are 

C(i_ij_i) = -dyz{xi, yj) ■ d,J{xi, yj)/ (2. ■ du ■ dv) 

C(i_ij) = (1 + dyz{xi,yjf)/du^ - dyZ{xi,yj) ■ d^z{xi,yj)/du 
C(i-ij+i) = dyz{xi,yj) ■ d^z{xi,yj)/{2 ■ du ■ dv) 
C(ij-i) = (1 + dxz{xi, yjf)/dv^ + da:z{xi, yj) ■ dxdyz{xi, yj)/dv 
-dyz{xi,yj) ■ dlz{xi,yj)/dv 
C(ij) = -2 ■ (1 + dyz{xi, yjf)/du^ - 2 ■ (1 + d^z{xi, yjf)/di 
C(ij+i) = (1 + d^z{xi, yjf)/dv^ - dxz{xi, yj) ■ dxdyz{xi, yj)/di 
+dyz{xi,yj) ■ dlz{xi,yj)/dv 
= dyz{xi, yj) ■ d^z{xi, yj)/{2 ■ du ■ dv) 
C(i+ij) = (1 + dyz{xi,yj)^)/du'^ + dyz{xi,yj) ■ d^z{xi,yj)/du 
-dyz{xi,yj) ■ dxdyz{xi,yj)/du 
C(i+ij+i) = -dyz{xi, yj) ■ d^z{xi, yj)/{2 ■ du ■ dv) (36) 

The right hand side is 

hd) = 2 • d^z{xi, yj) ■ d.^dyz{xi, yj) ■ dyz{xi, yj) - dlz{xi, yj) ■ (1 + dyz{xi, y^)^) 
- (1 + d^z{xi, yjf) ■ dp{xi, yj) (37) 

The linear equation 

Cdz'+\,j) = (38) 

for {N — 1) X {N — 1) length's vector corresponding to the points inside the 
boundary can be solved by using standard computer library. 

In the actual numerical calculation we multiply a reduction factor to the so- 
lution dz'^'^^ in each step to control the convergence. 
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4-2 Evaluation of the Area 



The standard expression [T] in the differential geometry for the area of a 
regular surface x('u, v) parameterized in terms of two scalar parameters u and 
V is 



area = / du / dv |xu x Xv|, (39) 



with 

c^x (9x 
x„ = — and Xt, = — . For x = x(x, y, z(x, y)), this becomes [1] 
ou ov 



area 



/ yr+^J+^dxdy, (40) 



where Q is the normal projection of the surface onto the xy plane. Accordingly, 
we calculated the area formed by the above mentioned discrete points as 



y/l + dxz{xi,yjY + dyz{xi,yjydu ■ dv. (41) 



This expression contains discretization errors. To estimate that, we discretized 
the bilinear interpolation for r = d = 1 in eq. (fT2|) as a 31 x 31 grid, and 
calculated the area obtained (of the discrete points) by this eq. fHTj) . This 
gave 1.2717 i.e. 0.7% underestimation of the exact value 1.280789 obtained by 
eq.dSni). 

In FigHl we show difference of the numerically calculated (N=40) minimal 
surface and the ruledi surface for r = d = 1. The corresponding difference of 
r = 2, d = 1 is shown in FigJSl 

That indicated that before reporting our 'numerical areas' we should compare 
different algorithms for calculating area out of a given set of points. Thus, we 
calculated the area by the sum of triangle 5*1 spanned by 

'"lij) = (O5 dv, z{ui, Vj) - z{ui, Vj-i)) and vf-^^ = {du, 0, z{ui, vj) - z{ui-i, vj)), 
and 5*2 spanned by ■y^j^) and v^ij) = {du,dv, z{ui,Vj) — z{ui-i,Vj-i)) 



E(l4.)x4,.)l + l4.)x4,.)l)/2 (42) 

(id) 
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Fig. 4. The difference of the numerical minimum and the ruled surface (r=l,d=l). 

The sum of triangles evaluated by the cross products is 1.281277037,i.e. 0.038% 
overestimation. 

The sum of triangles in the case of N=41 is 1.2811 i.e. 0.02% overestimation 
and in the case of N=21 is 1.2819 i.e. 0.09%. 

We also used a computer algebra system [T3] to find the two-dimensional 
interpolation surface working by fitting polynomial curves between successive 
data points followed by finding areas of the analytical interpolation surface 
X by an exact double integral of eq. (l39|l . The order 2 interpolation gave the 
above 1.280789195, the same up to 7 decimal places as the area without 

any discretization. 

Guided by this check, for areas formed by points we report both the areas cal- 
culated by triangulation as well by the interpolation-followed-by-the-double- 
integral; the numerical values strongly suggest these as better algorithms than 
the one used in eq. fHTll . 

The area of the ruled surface can be calculated analytically [lOj. In the Ap- 
pendix, we give formulae of the area of the ruledi surface and the ruled2 
surface. Numerically calculated area of the minimal surfaces( corresponding 
to the ruled2 surface) and analytically calculated area of the ruled surfaces 
for given r and d are compared in Table{T] The error bars are estimated from 
the convergence of the iteration. Numerical minimal surfaces corresponding 
to the ruledi are also slightly smaller than the analytical results. In the nu- 
merical calculation, approach to the absolute minimum is not guaranteed. In 
a variational calculation we could obtain slightly smaller area. 

An explicit analytical calculation of the minimal surface in R'^ is given in Ap- 
pendix 2. By constructing the conjugate minimal surface, K archer [TT] trans- 
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Fig. 5. The difference of the numerical minimum and the ruled surface (r=2,d=l). 





d 


numerical area 


ruled2 area 


ruledi area 


1, 


1 


1.2793(5) 


1.280789275 


1.280789275 


2, 


1 


2.3665(5) 


2.366974371 


1.861564196 


1, 


2 


3.1753(5) 


3.180414498 


4.316148066 


3, 


1 


3.4916(5) 


3.491711893 


2.595828045 


1, 


3 


5.9310(5) 


5.936348433 


9.325179471 


3, 


2 


7.2582(5) 


7.259880701 


6.208799631 


2, 


3 


8.5226(5) 


8.527786411 


10.22064879 



Table 1 

The numerical area (calculated using the order 2 interpolation) and the analytical 
area of the ruled2 surface for the hyperbolic paraboloid of given r and d. Analytical 
area of the ruledi surface is added for comparison. 

formed the plateau problem in into that in and showed that the global 
Weierstrass representation of triply periodic minimal surfaces is possible. We 
do not know whether the analytical calculation of the amount of the exact 
minimal surface area is possible through this method. 

We showed in the Appendix B that the exact minimal surface of Schwarz can 
be visualized. In order to evaluate the area, however, we need to interpolate the 
analytically obtained coordinates of the surface and perform numerial integra- 
tion. We leave this task as a future study. Accurate numerical evaluation of the 
amount of the area is important for physical application and the Chebyschev 
polynomial expansion is a practical method for performing this process since 
the area is parametrized as (x, z{x, y) instead of (a:(r, 9),y{r, 9), z{r, 9)). 
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A Appendix 1: Area of the ruled surface 



In this Appendix, we present the analytical formulae of the area of ruled 
surfaces [TOl. 



A.l Ruledi surface 



In the case of ruledi surface we define = (r, 0), = (— r, (i, 0), r'23 = 
(— r, 0, c?), ri4 = (r, 0, (i). The area is given by 

[ du j dv\[uri2 + (1 - M)r43] x [t>r23 + (1 - v)ri^]\ (A.l) 
Jo JO 



which becomes 



[\u C dvdJd^ + 2r2(l -2u + 2u'^) - Ar'^v{l - v) 
Jo Jo * 

= rf[Vrf2 + 2rV3 
-2rf^tan"^ 



-Ad^ + Adh-^ + 



+2d^tan~^ 

3d^ + P 
6r 

= d 

3rf2 + r2 



-Ad^ + Ad^r^ + 



/{12r') 
/{12r') 



log 



-r + Vd^ + 2r2 



6r 



log 



(A.2) 
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A. 2 Ruled2 surface 



The ruled2 surface is characterized by = {0,d,d), = (0,(i, — d), = 
(— r, 0, d), = (— r, 0, — rf). The area is given by 

du C dv\[urn + (1 - u)ri2] x [^23 + (1 - v)rii] \ (A.3) 
Jo JO 

which becomes 



du £ dvd^Jd^ + 2r2(l - 2u + 2u^) + _ 2vy 



dVd^ + 2r73 + — tan~H 
6 



^.dVd^ + 2r\ 



+ylog 



-d + + 



-r + Vd2 + 



(A.4) 



B Appendix 2: Visualization of the exact minimal surface 



In this Appendix, we construct conformal mapping from a complex uj plane to 
the skew quadrilateral of Schwarz, and visualize the surface using Mathematical^. 

The domain of the conformal mapping consists of an area bounded by four 

singula, points 6, c and i, where . = Zl±^, 6 = Zii^, , = iz^ 

v2 \2 y2 

1 - a/3 

and d = — — ^= — The Schwarz-Christoffel transformation corresponding to 
the four singular points would be expressed as 

R{uj) = f{uj)[{uj -a){uj- h){uj -c){uj- d)]-^'\ 

The Schwarz reflection principle implies, however, rotation of 180° about the 
boundary straight line is a symmetry of the mapping and the minimal surface 
area inside the boundary arc can be reflected to outside the boundary arc. 
Taking into account the presence of conjugate singular points, the actual Riu) 
is expressed as 

R{uj) = f{uj)[{uj -a){uj- b){uj -c){uj- d){uj - a'){uj - b'){uj - c'){uj - 

where a' = l/b,b' = 1/c, c' = 1/d and d' = 1/a. The position of the poles in 
the complex a; plane are given in Fig JB.li 
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We transform uj to ip, introduce a scaling parameter k and define 

2m 



R{zp) 



The coordinates of the minimal surface corresponding to the eq.([n]) scaled by 
K become 



2(1 + p2 



^ R f 



y 



—Im 



2(1 -p^ 



Vl + 14p4 + p^ 
p 4p 



— = Re , , , , 

K J VI + IV + 



(ip 
dp 
dp 



(B.l) 



The scaling parameter k, is defined at the end of the calculation. 



The boundary of the domain of the conformal mapping is bounded by four 
circles like 



UJ 



V2 



- <e< -. 
6 - - 3 



When 9 varies | ^ |, varies from ^^^^ y^ ^' ^ ^' 

The integral of x, z in the Weierstrass-Enneper representation given in sect. 2 
can be obtained by using the Mathematica, 



3 117 



- ^ [2pV/ - 4^3 + 7V/ + 4V/3 + 7 U , ^ _^ 



3Fi 



1115 



4'2'2'4' 7 + 4V3' _7 + 4V3, 



]/3Vp« + 14p4 + 1 



P 

(B.2) 



Im 111 








c: 




1 1 






















Re t 




o 


1 A 


d' 


















-0.5 -0-25 



Fig. B.l. Domains of the area in the 
complex UJ plane which are mapped to 
the Schwarz's minimal surface [1]. 



Fig. B.2. A ruled surface with the same 
boundary as that of Schwarz's minimal 
surface defined in [3]. 
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, _ + 4v/3 + 7)(/ - 4v/3 sinh'^ (-^) ||±|^) 

«^ 7(p8 + 14p4 + i)(7-4V3) 

(B.4) 

where Fi(a; 61, 62; c; x, is the Appell's 1st hypergeometric function, F{(j), m) 
is the eUiptic integral of the first kind. 



1 — i 

B. 1 The case of mapping inside the circle uj = ^ + V2e^^ 



1 — i 

The boundary of a circle whose center is at in the u plane {uj = re*") is 

given by 

1-i 



+ -\/2(cos 6* + z sin 6^) = r (cos a + i sin a) (B.5) 



V2 

We consider an area where a satisfies < a < 7r/2. The equation 

1 + 2cose 



cota = -, r = V3 + 2cose-2sin^ (B.6) 

-l + 2sin^' ^ ^ 

gives a solution of 9 as 

^, , _-, 1 cot^a + cot^a cot aV 3 — 2 cot a + 3 cot^ a , , 

%ota) = ±cos- i + eot^c ^ TT^^?^ 

(B.7) 

Here we choose the first sign + and the second sign - in the eq. (lB.7p . Then 
the r'^axi^^'^^) — 3 + 2 cos[6'(cot a)] — 2 sin[6'(cot a)] as a function of cot a 
behaves as Fig JB.3[ There is a branch point at cot a = 1, i.e. a = 7r/4. 



A mapping of a region < 6 < n/2 and < r < y/3 + 2 cos 6* — 2 sin 6' via 
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-Re 



2(1 + p2 



VI + 14p4 + pi 



:(ip 



i VI + 14p4 + p« ^ 



- = Re 



4p 



VI + 14p4 + p' 



-4p 



(B.8) 



is shown in Fig JB.4[ Due to the branch point near a = 7r/4, there appears 
numerical errors represented by thorns emanating from the saddle point. The 
blank area between the thorn going from the saddle point downwards and the 
left border of the minimal surface is due to numerical difficulties that inhibit 
simple extension of 6 and r to their boundaries. 



3 

U 2 
0- 




Fig. B.3. The maximal radius r^ax 
squared function of cota. In 

the calculation of the area, the region 
< cota is used. 




Fig. B.4. A piece of the Schwarz's min- 
imal surface. (The front right piece of 
the minimal surface of Fig. B.2.) 



B.2 The case of mapping inside the circle uo 



1 -i 



V2e'^(^+-/2) 



1 — i 

The boundary of the area of a circle whose center is at is given by 



1 - i 



+ \/2(cos(6' + 7r/2) + z sin(6' + 7r/2)) = r (cos a + i sin a) (B.9) 



The equation 



tana 



1 + 2cos(^ + 7r/2) _ l-2sine 
-1 + sin(^ + 7r/2) ~ -1 + 2cos^' 



V3 - 2cos^ 



2sin0 
(B.IO) 
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gives a solution of 6 as 



1 tan^a — tan^a tana 2 tana + 3 tan 



^(tana) = ±cos -[(1— tana ^ =F 



a. 



2 '^^ " " 1 + tan^ a 1 + tan^ a 

(B.ll) 

Here we choose both the first and the second sign + in the eq. flB.lip . The 
'"max (tan a) squared as a function of tana is shown in Fig |B.5[ 

A mapping of a region —n/2 < 9 < and < r < V3 — 2 cos 9 — 2 sin 9 via 



X 

K 

y 



-Re 



2(1 + p2 



Vl + 14p4 + p' 



-4p 



K J VI + IV + 



2; 



4p 



Vi + IV + 



dp 



(B.12) 



is shown in Fig JB.6t which gives a part of the minimal surface. There appears 
a missing region in upper left corner due to the singularity of tana near 
a = 90°. The small missing region is mapped in a place rotated by 90° and 
z < —2 X 0.47196 region as shown in Fig JB.7[ 

To evaluate the area of the minimal surface, we use data of the right half of 
the triangle of Fig JB.7[ 




Fig. B.5. The maximal radius r^ax 
squared function of tana. In 

the calulation of the area, the region 
tan Q < is used. 




Fig. B.6. A piece of the Schwarz's min- 
imal surface (The front left piece of the 
minimal surface of Fig. B.2. ) 



A combination of the Figs lB.4l and IB. 61 are shown in Fig JB.81 The scale pa- 
rameter K is defined by the height of the coordinate of the edge of the 
tetrahedron 
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Fig. B.7. The bird's eye view of the 
minimal surface of Fig lB.61 The miss- 
ing corner region in that figure appears 
in the 90° rotated region. 

r(v^-l)i/V2 4p 

-dp 



Jo 



2i (2 + Vs) sinh"^ (? 




Fig. B.8. A combination of the two 
pieces of the Schwarz's minimal surface. 



4^3) I97 + 56V3) = -0.47196 (B.13) 



The Fig lB.81 indicates that the actual height of the z- 
this value and since this height should be — ^, k, 



0.37456. 



coordinate is twice of 
^/(2 X 0.47196) = 



We observe that the scale factor given in Ref. [1] does not agree with ours. 
To the best of our knowledge, it is the first explicit calculation of the exact 
minimal surface whose bounding contour is the skew quadrilateral. 
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